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摘 要 : 及 时 获取 凌 汛 期 河 冰 和 水 体 的 空间 分 布 特征 ,对 于 预测 冰凌 灾害 、 提 高 防凌 信息 化 管理 水 


平 有 重要 意义 。 遥 感 技术 是 当前 获取 河 冰 和 水 体 空间 分 布 的 最 主要 手段 之 一 。 但 是 ,黄河 水 体 有 
量 悬 浮 泥 沙 ,这 给 基于 蜗 感 技术 的 高 精度 冰 - 水 分 类 带 来 了 挑战 。 以 黄河 内 蒙古 段 为 例 ,基于 


Landsat 8 OLI 遥 感 影像 数据 ,在 利用 归 一 化 积 雪 指 数 (NDSI) 及 河道 矢量 数据 排除 无 关 地 物 的 基础 
上 ,对 比 了 近 红 外 波段 反射 率 值 \ 归 一 化 差异 水 体 指数 (NDWI)、 归 一 化 积 雪 指 数 (NDSI) 改进 的 归 
一 化 积 雪 指数 (MNDSI) 以 及 归 一 化 差异 未 封冻 水 体 指 数 (NDUWI) 在 黄河 内 蒙古 段 典 型 河道 河 冰 、 
水 体 分 类 中 的 表现 ,计算 各 指标 总 体 分 类 精度 及 Kappa 系 数 并 进行 国 值 稳 定性 分 析 。 结 果 表 明 : 在 
利用 NDSI 和 高 清 历 史 影像 排除 河道 外 无 关 地 物 后 ,NDUWI 在 各 子 段 影像 中 的 总 体 分 类 精度 和 
Kappa 系 数 均 达到 90.00% 及 0.90 以 上 ,其 河 冰 \ 水 体 最 优 区 分 国 值 大 体 分 布 于 冰 值 中 值 附近 。 研 究 
结果 可 为 凌 汛 期 黄河 冰凌 监测 方法 的 选取 以 及 冰 上 爆破 位 置 的 拟定 提供 依据 。 

X 键 词 : 河 冰 ; 遥感 指数 ， 阅 值 稳定 性 ; 不 确定 性 ; Landsat8; 黄河 内 蒙古 段 
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冰凌 由 河流 中 的 河水 冻结 .堆积 而 产生 ,在 河 
道中 堆积 易 引 发 冰 塞 ` 冰 坝 等 现象 , 轻 则 河 冰 堆 积 
EE, HUE ESE .省 堤 决口 。 黄 河内 蒙古 段 河 
道 形 态 多 样 , 河 势 复 杂 , 易 产生 冰凌 灾害 ”2 。 因 此 ， 
准确 提取 河 冰 ,水体 信息 对 研究 河 冰 生 消 过 程 . 河 
冰 空 间 变 化 以 及 冰凌 洪水 的 研判 至 关 重 要 。 传 统 
野外 观测 局 限于 “点 对 点 ” 式 的 局 部 观测 ,无 法 获得 
完整 的 河 冰 信 息 , 冰 情 观测 受到 限制 。 卫 星 遥 感 则 
具有 观测 范围 大 .时效 性 强 .数据 内 容 客观 等 优点 ， 
近 些 年 在 河 冰 研究 领域 已 被 广泛 应 用 3? 。 

日 前 ,利用 遥感 技术 提取 河 冰 、 水 体 信息 主要 
以 目 视 解 译 | 0 Ez c 58] S BE. . 归 一 化 积 雪 指数 
(NDST) 、 归 一 化 差异 水 体 指 数 (NDWI) 以 及 改进 的 
归 一 化 积 雪 指数 (MNDSI) 与 单 波段 靖 值 相 结 合 的 方 
法 为 主 。 杨 中 华 等 基于 “四 星 三 源 ” 模 式 , 在 国 
内 最 早 将 遥感 技术 应 用 于 黄河 冰 情 监测 当中 ,限于 
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当时 卫星 影像 空间 分 辨 率 , 作 者 仅 进 行 大 尺度 宏观 
观测 ,并 未 进行 局 部 河 段 微观 分 析 ; 赵 水 霞 等 基于 
Landsat 8 彩色 合成 影像 对 黄河 什 四 份子 弯 道 河 冰 
过 程 及 分 布 进行 目 视 解 译 ,但 作者 并 未 进行 提取 精 
度 验证 ;Chu 等 "基于 MODIS 数 据 近 红外 波段 监测 
了 加 拿 大 奴隶 河 河 冰 过 程 , 作 者 并 未 确定 河 冰 、 水 
体 在 近 红 外 波段 的 反射 率 值 , 而 是 假设 当 近 红外 波 
段 的 反射 率 最 大 时 河道 内 均 为 河 冰 ,而 反射 率 最 小 
时 河道 内 均 为 水 体 ;Chaouch 等 “利用 MODIS 可 见 
光 及 近 红 外 波段 的 国 值 构建 决策 树 ,提取 河 冰 信息 
并 获取 河 冰 范 围 ;Dozier” 于 1989 年 提出 的 用 于 监 
测 积 雪 变 化 的 雪 信息 提取 NDSI 指 数 ,该 指数 对 积 雪 
与 土壤 的 区 分 度 较 高 ,目前 也 被 应 用 于 大 尺度 冰 
川 、 湖 冰 等 领域 * 呈 ;Li 等 外 利用 NDSI 对 1999 一 
2018 年 青藏 高 原 八 宝 河 河 冰 分 布 进行 监测 ,而 八 宝 
河流 域 冬季 河水 完全 封冻 并 不 存在 清 沟 ( 狭 长 型 不 
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结 冰 的 裸露 水 面 ) ,提取 任务 主要 以 区 分 河 冰 与 土 
壤 像 元 为 主 , NDSI 分 类 精度 令 人 满意 ; 件 献 友 等 
根据 黄河 水 体 的 光谱 特征 ,以 蓝 色 及 绿色 波段 构建 
归 一 化 差异 未 封冻 水 体 指数 (NDUWI) ,在 利用 
NDSI 排 除 土 壤 、 树 木 等 其 他 地 物 像 元 的 基础 上 , 进 
而 利用 NDUWI 进 行 河 冰 水 体 提 取 , 但 以 水 体 作为 切 
人 和信 点 的 提取 方法 的 精度 有 待 进一步 验证 ; 勾 鹏 等 “发 
现 近 红外 波段 与 短波 红外 波段 构建 的 归 一 化 指数 
直方 图 中 两 峰值 间距 大 于 NDSI 直 方 图 中 的 两 峰 间 
距 , 有 利于 区 分 河 冰 水体 ,由 此 建立 改进 的 归 一 化 
差异 积 雪 指数 (MNDSI) 对 纳木错 湖 冰 进行 识别 , 取 
得 很 好 的 效果 。 综 合 来 看 ,对 于 黄河 河 冰 、 水 体 的 
监测 ,当前 遥感 技术 和 卫星 影像 的 时 空 分 辨 率 已 基 
本 满足 冰期 黄河 冰凌 监测 的 能 力 ,但 单一 使 用 遥感 
旨 数 模型 存在 难以 区 分 河 冰 \ 水 体 与 土壤 的 问题 ， 
河道 边界 难以 提取 ;同时 ,各 指数 模型 对 于 黄河 河 
冰 、 水 体 监 测 的 精度 还 有 待 验 证 。 
黄河 内 蒙古 段 由 于 其 自身 河道 水 力 条 件 的 特 
殊 性 , 凌 汛 期 易 发 生 卡 冰 结 坝 现象 ,增加 “ 武 开 河 ” 
的 风险 ”。 鉴 于 此 ,本 文选 取 黄河 内 蒙古 段 5 个 试 
验 河 段 为 研究 对 象 , 基 于 Landsat 8 OLI 遥 感 影像 数 
Jis ,在 利用 Google Earth 历史 高 清 影像 结合 NDSI 指 
数 模型 排除 河道 外 无 关 地 物 的 基础 上 ,使 用 近 红 外 
波段 反射 率 值 、 NDWI、NDSI、MNDSI LA & NDUWI, 
利用 各 指数 闵 值 对 河 冰 水体 进行 分 类 ,进而 计算 
监测 精度 ,对 立 值 的 稳定 性 及 方法 的 适用 性 进行 评 
价 , 以 期 为 遥感 技术 监测 河流 冰 情 及 冰凌 洪水 的 研 
判 提 供 参 考 。 


1 研究 区 数据 与 方法 


1.1 研究 区 概况 

研究 区 始 于 黄河 内 蒙古 段 海 勃 湾 水 库 ,终于 万 
家 寨 水 库 库 尾 区 ,全 长 约 为 650 km, E [fc 3" JL" 5e 
形 分 布 呈 (图 1)。 分 别 在 河 段 上 中、 下 游 选 取 试 验 
河 段 , 河 段 类 型 包括 过 渡 型 ,游荡 型 及 弯曲 型 ;过 渡 
型 河 段 浅滩 多 日 河道 摆动 幅度 大 ;游荡 型 河 段 淤积 
PR AA RA VIBE (EES MA, SUK 
Py Se Pos 798 £355 TERRE ok TREWEDK ` 清 沟 及 土壤 , 河 冰 
主要 以 平滑 冲击 岸 冰 为 主要 类 型 。 河 段 每 年 12 月 
初 开始 流 凌 ,至 翌年 3 月 未 河 段 冰凌 全 部 消融 , 受 来 
水 条 件 .河道 特性 等 因素 影响 ,水 体 在 河床 内 及 河 


— 图 例 
0 30km 水 流 流 向 全 水 文 站 
一 黄河 


图 1 研究 区 概况 
Fig. 1 Overview ofthe study area 


滩 上 结 冰 ;在 流速 较 大 的 急流 ,浅滩 处 或 排水 口 处 
常 出 现 狭长 型 不 结 冰 的 裸露 水 面 ( 清 沟 )。 
1.0 数据 来 源 与 处 理 
1.2.1 Landsat 8:& & X 4& — 河 冰 信息 提取 基于 
Landsat 8 OLI 影 像 ( 后 文 简称 L8 影 像 ) 数 据 , 数 据 获 
取 自 美国 地 质 调查 局 官网 (United States Geological 
Survey ,简称 USGS ,https://earthexplorer.usgs.gov/) ,所 
用 影像 均 为 LIT 级 产品 ,该 级 别 影像 是 使 用 地 面 控 
制 点 和 数字 高 程 模型 数据 进行 精确 校正 的 数据 产 
品 ,误差 小 于 0.4 个 像 元 。 试 验 所 用 使 用 的 影像 为 
无 云 或 河道 及 河道 附近 无 云 的 影像 数据 。 影 像 信 
息 见 表 1。 
1.2.2 验证 数据 试验 使 用 Google Earth 历史 高 清 
影像 作为 参考 影像 ,选取 与 L8 影 像 拍摄 日 期 同日 或 
相近 日 期 的 历史 高 清 影 像 ,由 于 拍摄 日 期 处 在 河 冰 
稳定 封冻 的 1 一 2 月 ,L8 影 像 与 参考 影像 成 像 时 间 差 
内 的 河 冰 履 盖 基本 不 变 。 基 于 参考 影像 目 视 解 译 
出 的 各 类 主要 地 物 矢 量 数据 将 用 于 分 类 精度 评 
价 。 因 参考 影像 的 空间 分 辩 率 在 亚 米 级 , 故 目 视 解 
译 的 误差 将 忽略 不 计 。 
1.2.3 数据 预 处 理 L8 影 像 的 预 处 理 过 程 为 辐射 定 
标 及 大 气 校正 。 预 处 理 过 程 利用 集成 在 ENVI 5.5 
版 本 中 的 辐射 定 标 及 FLAASH (Fast line-of-sight 
atmospheric analysis of spectral hypercubes ) 大 气 校正 
工具 获得 地 表 反 射 率 数 据 。 其 中 ,FLAASH 中 大 气 
模型 参数 根据 影像 拍摄 日 上 午 的 大 气 柱 水 汽 量 确 
定 ,该 数据 获取 自 MODIS 标 准 大 气 产 品 中 的 MOD05 
数据 集 ( 下 载 地 址 :http:Wladsweb.nascom.nasa.gov) 。 
大 气 校正 参数 见 表 1。 

所 有 L8 影 像 均 经 过 堤防 矢量 数据 进行 裁剪 ,只 
保留 两 岸 堤防 间 的 部 分 。 首 先 基 于 Google Earth Jj 
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表 1 Landsat 8 OLI 遥 感 影 像 信息 及 大 气 校正 参数 


Tab.1 Landsat 8 OLI remote sensing image information and atmospheric correction parameters 


Landsat 8 OLI 日 期 平均 大 气 柱 


Landsat 8 OLI 日 期 FHK a 


影像 条 带 号 。 (年 -月 -日 KRipo 。 混合 层 高 度 em | 影像 条 带 号 (年 -月 -日 ) KRipo PARRE 

129.032 2014-01-17 0.29 1.42 129.032 2017-02-26 0.22 2.19 
128.032 2016-02-01 0.19 1.69 127.032 2018-02-15 0.12 2.23 
128.032 2018-02-22 0.10 1.20 128.032 2018-02-06 0.19 2.20 
130.032 2020-02-10 0.13 1.43 129.032 2018-02-13 0.17 2.23 
128.032 2020-02-28 0.22 0.96 127.032 2019-01-17 0.10 2.19 
127.032 2015-02-23 0.17 2.10 128.032 2019-01-24 0.26 2.13 
128.032 2015-02-14 0.12 2.12 129.032 2019-02-16 0.28 2.22 
129.032 2015-02-21 0.23 2.13 127.032 2020-02-21 0.22 2.13 
127.032 2016-02-10 0.21 2.22 128.032 2020-02-12 0.14 2.13 
128.032 2016-02-01 0.22 2.13 129.032 2020-02-19 0.22 2.22 
129.032 2016-02-08 0.15 2.19 127.032 2021-01-22 0.29 2.22 
127.032 2017-02-28 0.11 2.14 128.032 2021-01-29 0.29 2.19 
128.032 2017-02-19 0.13 3.65 129.032 2021-02-21 0.23 2.13 

史 高 清 影像 目 视 解 译 绘制 的 河道 矢量 边界 进行 裁 人 

DY ,排除 大 部 分 河道 外 无 关 地 物 , 其 次 利用 NDSI 对 i 

河岸 边缘 土壤 及 其 他 无 关 地 物 进行 二 次 排除 。 大 气 校正 

实验 区 位 于 两 岸 堤防 内 的 河道 范围 ,相关 法 律 规 定 TY 

堤防 内 部 不 得 设立 建筑 ,道路 等 设施 且 河 滩地 不 得 Eun 

占用 , 故 经 堤防 矢量 裁剪 后 的 影像 数据 可 有 效 排除 


附近 民居 .道路 .工厂 等 其 他 无 关 地 物 ; 又 因 河道 内 
树木 极 少 且 处 于 冬季 , 故 假 定 河道 内 仅 有 河 冰 水 
体 .土壤 3 种 地 物 。 而 NDST 对 土壤 和 冰雪 的 分 类 精 
度 已 有 前 人 证 实在 此 不 再 著述 。 为 尽 可 能 保留 全 
部 河 冰 .水体 像 元 , 故 以 0 作为 土壤 与 河 冰 、 水 体 的 
分 类 阔 值 ,影像 预 处 理 及 分 类 流程 见 图 2 ,流程 化 处 
理 后 待 提取 像 元 仅 为 河 冰 与 水 体 ,几乎 无 其 他 无 关 
地 物 。 
1.3 研究 方法 
1.3.1 5 种 适 感 指数 模型 ”选取 5 种 壳 感 指数 模型 
对 黄河 内 蒙古 段 5 个 子 段 的 河 冰 水体 进行 分 类 ,分 
别 为 : 

(1) 近 红 外 波段 反射 率 阔 值 法 "。 水 体 在 近 红 
外 波段 反射 率 值 低 ,而 冰雪 的 反射 率 较 高 , 故 适 用 
于 冰雪 与 水 体 的 分 类 。 分 类 算法 如 下 : 

; AMA, pn/10000 < threshold 
pee Ir p.,,/10000 > al Wd 

式 中 : pyr 为 近 红 外 波段 反射 率 值 ; pixel HERL 
像 像 元 点 ;threshold 为 浆 值 ,下 同 。 为 节省 储存 空 
HJ , USGS 发 布 L8 影 像 数 据 的 像 元 反射 率 值 均 放 大 


OU -zecuz 


l 
—3 


总 体 分 类 精度 、 
Kappa 系 数 


注 :NDWI 表 示 归 一 化 差异 水 体 指数 ;NDSI 表 示 归 一 化 差异 积 雪 指 

数 ;MNDSI 表 示 改 进 的 归 一 化 差异 积 雪 指数 ;NDUWI 表 示 归 一 化 
未 封冻 水 体 指数 。 下 同 。 

图 2 影像 预 处 理 及 分 类 实验 流程 图 


Fig.2 Flowchart of image preprocessing and classification 


experiment 


10000 倍 ,计算 时 除 以 10000 将 取 值 范围 限制 在 0 到 
1 之 间 以 便于 与 其 他 指数 值 进行 比较 。 

(2) 归 一 化 差异 水 体 指 数 (NDWI) 阅 值 法 。 
NDWI 利 用 绿色 波段 与 近 红 外 波段 构建 归 一 化 差异 
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指数 ,是 由 MecFeeters2 提 出 的 目前 最 为 常用 的 陆 表 
水 体 指 数 。 因 河道 内 并 不 涉及 建筑 用 地 ,NDWI 的 
缺点 不 会 影响 其 对 河 冰 、 水 体 提 取 。 分 类 算法 
如 下 : 
az JAR (Pore T Pi poss +P yn) > threshold 
pes oi (Puna Dads bis ud 
(2) 
RP: pu. 为 绿色 波段 反射 率 值 ,下 同 。 
(3) 归 一 化 积 雪 指数 (NDSI) 阔 值 法 NDSI 
是 由 绿色 波段 与 短波 红外 波段 构建 的 归 一 化 差异 
指数 ,此 外 其 与 MNDWI 指 数 公 式 一 致 。MNDWI 由 
PAK HE , 因 水 体 在 SWHR 波段 的 反射 率 值 更 
低 , 较 NDWI 可 有 效 区 别 水 体 与 其 他 地 物 。 分 类 算 
法 如 下 : 
，， | 水 体 , (p6 Pas (Poss t Pss) > threshold 
pe ( (po Peon Gea + pswm) < e 
(3) 


式 中 : oss 为 短波 红外 波段 反射 率 值 ,下 同 。 

(4) 改进 的 归 一 化 差异 积 雪 指 数 (MNDSI)P 。 
MNDSI 利 用 近 红 外 及 短波 红外 波段 构建 归 一 化 差 
异 指数 ,而 其 计算 方法 与 Ouma 等 中 提出 的 系列 水 
体 指 数 中 的 NDWE 计 算 方 法 一 致 ,该 方法 对 水 体 边 
缘 提取 效果 较 好 。 分 类 算法 如 下 : 

，， | 水 体 , (Pswm Pus (Pss +Pun) > threshold 
piel - [ru (Pss — Pra) Pew +P xm) > ene 
(4) 

(5) 归 一 化 差异 未 封冻 水 体 指 数 (NDUWD 浆 值 
法 号。 件 献 友 等 扩 提出 的 NDUWI 利 用 蓝 绿 波 段 构 
建 归 一 化 差异 指数 ,并 利用 该 指数 对 黄河 内 蒙古 段 
河 冰 进 行 监测 。 

| PRIE, (Poren Prine) (ss + P) > threshold 
puer m (Dcus “Pre Poss t Pn) € bd 
(5) 


式 中 : py, 为 蓝 色 波段 反射 率 值 。 

1.3.2 精度 评价 为 比较 各 算法 分 类 精度 ,利用 总 
体 分 类 精度 (po, 指 对 每 一 个 随机 样本 ,所 分 类 的 结 
果 与 检验 数据 类 型 相 一 致 的 概率 ,单位 为 %) 和 Kap- 
pa 系数 (x, 基 于 混淆 矩阵 衡量 河 冰 、 水 体 分 类 精度 
的 指标 , 介 于 [-1, 1]) 评 价 各 分 类 方法 的 提取 精度 
(p, 为 中 间 变 量 ,表示 偶然 一 致 性 ,单位 为 %) ,其 计 
算 过 程 见 式 (6)~(8)。 根 据 各 方法 的 提取 结果 列 出 
河 冰 \ 水 体 混淆 矩阵 ( 表 2) ,并 分 别 计算 各 指数 分 类 


表 2 河 冰 \ 水 体 混 淆 矩阵 


Tab.2 Confusion matrix of river ice and water body 


河 冰 真 实 值 水 体 真实 值 
河 冰 预测 值 Gu an 
水 体 预测 值 az an 


HJ po ek. Landis E f 120.81 时 分 类 几乎 完全 
吻合 ,但 为 达到 精确 提取 的 目的 ,选取 x>0.9 作 为 分 
类 达到 几乎 完全 吻合 的 标准 。 
Po = + ala, + Ay, + Ay, * Ay (6) 
P.= {lan + a4) X(a,, a4) + (as, + a5] x 
(an + ay)Ma,, + a4, * a4, tan) 


K 7 (po 7p.) (1 - p.) (8) 
2 ”结果 与 分 析 


2.1 河 冰 、 水 体 分 类 结果 

分 别 采用 5 种 遥感 指数 模型 对 黄河 内 蒙古 段 河 
冰 \ 水 体 信息 进行 提取 ,限于 篇 幅 , 仅 展示 2016 年 和 
2018 年 2 期 影像 提取 结果 (图 3、 图 4)。 总 体 来 看 ， 
基于 NDSI 结 合 Google Earth 历史 高 清 影像 提取 到 的 
河道 边界 排除 了 土壤 及 河岸 滩地 上 冰雪 等 地 物 , 降 
低 了 无 关 地 物 对 于 后 续 总 体 分 类 精度 及 Kappa 系数 
计算 精度 的 影响 ,从 分 类 效果 来 看 ,5 种 方法 分 类 效 
果 由 高 到 低 依 次 为 : 归 一 化 未 封冻 水 体 指 数 、 归 一 
化 差异 水 体 指 数 、 近 红外 波段 反射 率 值 改进 的 归 
一 化 差异 积 雪 指 数 . 归 一 化 差异 积 雪 指数 。 在 2018 

影像 中 ,NDSI 较 其 余 4 种 方法 在 河 冰 、 水 体 信息 

提取 上 欠 佳 (图 3), 误 分 现象 严重 ;NDWI 与 NDUWI 
在 遥感 影像 中 分 类 效果 最 佳 ,其 中 NDUWI 总 体 分 
类 精度 和 Kappa 系 数 均 达 到 95.00% 和 0.90 以 上 。 

从 5 种 遥感 指数 模型 所 对 应 的 最 大 (最 小 ) 国 值 
来 看 ( 表 3) , 河 冰 \ 水 体 分 类 结果 与 阀 值 大 小 并 无 直 
接 关 系 ,NDSI 阅 值 范围 最 大 ,但 其 分 类 结果 在 5 个 
试验 区 并 不 稳定 ,NDWI 与 NDUWI 阔 值 范围 适中 ， 
但 分 类 效果 在 5 个 试验 区 内 分 类 结果 较 好 。 进 一 步 
分 析 发 现 ,NDUWI 在 5 期 遥感 影像 中 的 河 冰 、 水 体 
最 优 区 分 阔 值 大 体 分 布 于 阔 值 中 值 附 近 , 如 在 2016 
年 遥感 影像 中 NDUWI 的 浆 值 为 (0~0.46) ,最 优 效 值 
为 0.21。 
2.2 阅 值 稳定 性 分 析 
2.2.1 敏感 性 分 析 图 5 给 出 了 试验 区 5 期 遥感 影 
像 河 冰 水 体 分 类 精度 对 于 阔 值 变化 的 响应 。 从 各 


(7) 
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(b EUDUDPID EXE oa 


SRH 水 体 Yalu 


注 :m 指 对 每 一 个 随机 样本 ,所 分 类 的 结果 与 检验 数据 类 型 相 一 致 的 概率 ;k 是 Kappa 系 数 。 下 同 。 
图 3 2016 年 2 月 1 日 5 种 遥感 指数 模型 河 冰 水 体 分 类 结 


Fig. 3 Classification results of ice-water in five remote sensing models on February 1, 2016 


期 影像 中 各 方法 的 分 类 精度 来 看 ( 表 4) ,在 2020 年 
的 2 期 影像 中 近 红 外 波段 反射 率 阔 值 法 在 相 较 其 他 
方法 展现 出 最 高 的 提取 精度 ,m 及 x 分 别 为 99.18% 、 
0.99 及 99.98% .0.99;2014 .2016 年 及 2020 年 影像 中 
NDUWI 阅 值 法 相 较 其 他 方法 展现 出 最 高 的 分 类 精 
FE , po Be HI Ay 99.38% , 98.90% , 99.8096 及 0.98 , 
0.97、0.99。 从 各 方法 在 各 期 影像 中 的 分 类 表现 来 
看 , 近 红 外 波段 反射 率 阔 值 法 在 2018 年 1 期 和 2020 
年 2 期 共 3 期 影像 中 po 及 k 均 达到 并 超过 了 90.00% 
及 0.90,2014 年 和 2016 年 影像 中 po 未 达到 90.00% , K 

ARAL 0.7 ; Eje £L CE sc] SS LAL, NDWI 
在 前 2 期 影像 中 分 类 效果 不 理想 ,在 后 3 期 影像 中 p。 
及 k 均 达到 并 超过 90.00% 及 0.90; NDSI 在 2014 及 
2016 年 2 期 影像 中 ps 及 xk 均 达到 并 超过 90.00% 及 
0.90, 而 在 后 3 期 影像 中 «出 现 为 0 的 情况 ,说 明 误 
分 . 漏 分 现象 严重 ; MNDSI 在 5 期 影像 中 po 及 kk 总体 


Ti ETE 85.00% x 0.80 ÆA ; NDUWI Bü (EA TE 5 HH 
感 影像 中 的 po 及 k 均 达到 并 超过 了 90.00% 及 0.90。 
进一步 分 析 发 现 ( 图 5), 当 采用 适当 阔 值 时 ,对 
于 上 述 5 种 遥感 指数 模型 ,所 能 达到 的 最 高 精度 从 
低 到 高 排序 依次 为 :NDSI、MNDSI、 近 红外 波段 反射 
率 阐 值 法 .NDWI.NDUWI。 此 外 ,从 分 类 精度 对 冰 
值 变 化 的 响应 机 制 来 看 , 近 红 外 波段 反射 率 阔 值 法 
变化 幅度 最 大 , 即 在 分 类 时 只 要 姜 值 有 微小 的 变 
动 ,就 会 造成 提取 结果 的 显著 变化 ,总体 精度 差 值 
为 50.10%; 相 比较 而 言 ,NDWI 阔 值 变化 对 分 类 精度 
敏感 程度 较 低 ,选取 恰当 阔 值 , 便 可 得 到 较 好 的 精 
度 ;对 于 其 他 3 种 指数 模型 , 闷 值 对 分 类 精度 影响 的 
敏感 程度 由 低 到 高 分 别 为 :NDSI、MNDSI、NDUWI。 
总 体 来 看 ,NDUWI 对 阔 值 变化 的 敏感 程度 较 小 ,总 
1k ps Bert MNDSI my , fr 3 “4 ELAN po Rew E Bem 
100% 和 1, 分 类 精度 高 。 
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图 4 2018 年 2 月 22 日 5 种 遥感 指数 模型 河 冰 水 体 分 类 结果 


Fig.4 Classification results of ice-water in five remote sensing models on February 22, 2018 


RI 5 种 遥感 指数 模型 对 应 最 大 (最 小 ) 阅 值 范 围 


Tab.3 Maximum (minimum) threshold range of the five remote sensing models 


BY (Eye Hl 近 红外 波段 反射 率 值 NDWI NDUWI NDSI MNDSI 
最 小 0~0.37 0~0.37 0~0.22 0~0.86 -0.85~-0.05 
最 大 0~0.62 0~0.99 0~0.78 0~1.00 -0.98-0 


雪 指 数 。 下 同 。 


综合 试验 区 5 期 提取 结果 发 现 ,L8 影 像 不 同 指 
数 模型 在 同一 研究 区 内 阔 值 对 河 冰 水 体 提取 精度 
的 影响 程度 不 同 ;相同 指数 模型 在 不 同 研究 区 内 所 
能 达到 的 提取 精度 不 同 , 闵 值 对 河 冰 .水体 提取 精 
度 影响 程度 也 不 同 ; 例 如 ,在 2016 年 遥感 影像 中 , 当 
选用 一 定 阀 值 时 ,NDUWI 精 度 最 高 ,MNDSI 次 之 ;而 
在 2018 年 影像 中 ,NDUWI 在 一 定 阐 值 内 精度 最 高 ， 


注 :NDWI 表 示 归 一 化 差异 水 体 指数 ;NDUWI 表 示 归 一 化 未 封冻 水 体 指 数 ;NDSI 表 示 归 一 化 差异 积 雪 指数 ;MNDSI 表 示 改 进 的 归 一 化 差异 积 


AL NR E BI SS B EK o 

2.2.2 不 确定 性 分 析 基于 遥感 影像 数据 对 黄河 内 
蒙古 段 冰 水 信息 提取 中 ,由 于 主观 和 客观 因素 的 存 
在 ,往往 影响 到 最 后 的 提取 精度 。 例 如 ,NDSI 在 
2014 年 和 2016 年 2 期 影像 中 分 类 精度 较 高 ,而 在 后 
3 期 影像 的 分 类 上 差强人意 。 本 研究 的 不 确定 性 来 
源 于 以 下 几 个 方面 :遥感 影像 的 选取 人工 目 视 解 
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~ MNDSI  — NDUWI 
12 T (b) 2014 年 1 月 17 日 Kappa 系 数 
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图 5 IEE po Alla Xt EJIE ÉS fs Jz 
Fig. 5 Response of po and « in the test area on the threshold 
表 4 5 种 遥感 指数 模型 对 应 最 高 分 类 精度 
Tab.4 Highest classification accuracy of the remote sensing index models 
分 类 精度 近 红 外 波段 反射 率 值 NDWI NDUWI NDSI MNDSI 
pol% 99.12 98.25 99.89 98.34 99.2] 
K 0.98 0.98 0.99 0.9 0.98 


注 :po 指 对 每 一 个 随机 样本 ,所 分 类 的 结果 与 检验 数据 类 型 相 一 致 的 概率 ;x 是 Kappa 系数。 
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FHEWY 


译 经 验 及 研究 区 环境 变化 。 冬 季 , 由 于 河道 上 空 
筋 云 及 手 的 存在 ,光学 数据 在 冰 情 监测 中 可 用 数 
据 量 有 限 , 每 个 试验 区 仅 能 王选 1~2 景 可 用 于 试验 
的 影像 汪 。 本 研究 所 采用 的 L8 影 像 均 为 无 云 或 河 
道 及 河道 附近 无 云 的 影像 数据 ,尽量 保证 试验 数据 
的 质量 及 可 靠 性 。 受 客观 因素 影像 ,目前 尚未 引入 
其 他 卫星 影像 作为 补充 数据 源 。 人 工 目 视 解 译 的 
误差 主要 来 自 混 合 像 元 ,利用 阀 值 模型 计算 的 栅 格 
数据 一 般 包 含 多 种 地 物 信息 , 除 河道 内 部 少量 平 请 
冰 外 %( 即 静态 岸 冰 , 指 平均 流速 低 于 河道 临界 流 
速 的 水 力 条 件 下 ,在 河道 内 同时 沿 纵 向 和 横向 发 展 
的 表面 较为 光滑 的 河 冰 ), 少 有 纯净 像 元 ,特别 是 河 
道 边 界河 冰 水 体积 雪 以 及 土壤 之 间 界 限 难 于 准 
确 拟 定 , 在 人 为 判别 时 易 影 响 分 类 精度 ”。 在 变化 
环境 方面 ,大气 中 水 汽 CARR RA .河道 清 沟 内 
水 体 的 含 沙 量 、 河 岸 边界 土壤 湿度 等 造成 “ 同 谱 异 
质 、 同 质 异 普 ” 等 现象 ,对 边界 提取 、 冰 水 情 监 测 造 

i 影响”。 为 了 最 大 程度 上 降低 河道 边界 提取 中 
的 误 判 ,本 研究 设计 了 NDSI 结 合 Google Earth 历史 
高 清 影像 对 边界 提取 的 流程 化 提取 操作 ,使 得 土壤 
等 无 关 地 物 最 大 程度 得 以 排除 。 


ua 


3 讨论 


黄河 内 蒙古 段 全 长 约 为 650 km, 东西 跨度 较 
大 ,整体 呈 “ 几 ”字形 分 布 ,野外 冰 情 受到 限制 ,基于 
逐 感 影像 在 冰 情 监测 方面 效果 较 她 。 同 时 ,黄河 
冰 情 信息 监测 受 限 于 数据 的 时 相 分 辨 率 以 及 天 气 
条 件 , 因 此 ,党 在 研究 中 出 现 某 一 稳 封 期 数据 缺失 
的 情况 “。 与 传统 的 大 尺度 海 冰 ( 湖 冰 ) 冰 情 监测 不 
同 ,由 于 黄河 自身 的 泥 沙特 性 以 及 河道 洲 积 昌 延 情 
况 ,本 文 研究 基于 L8 影 像 ,综合 Coogle Earth 历史 高 
清 影 像 和 NDSI 首 先 排 除了 河道 外 士 壤 、 建 筑 及 植被 
等 无 关 地 物 , 仅 保留 冰 水 信息 以 保证 被 研究 像 元 的 
真实 性 ,解决 了 河道 边界 河 冰 水体 信息 难以 识别 
的 问题 3。 

本 研究 发 现 ,遥感 影像 选取 .河道 上 空气 象 条 
件 以 及 河流 水 质 条 件 是 影响 冰 情 监测 提取 精度 的 
重要 因素 。 在 数据 选取 方面 ,由 于 冬季 河道 上 空 
筋 云 及 手 的 存在 ,不 同 指数 分 类 阔 值 难以 拟定 , 存 
在 较 大 浆 值 波动 情况 , 故 应 选用 稳 封 期 ( 即 河 冰 不 
再 治 河道 横向 发 展 ) 无 云 . 雪 且 气 象 条 件 恨 好 时 期 


的 卫星 影像 站。 在 冰 水 信息 提取 实验 设计 方面 ， 
NDSI 哩 不 能 精确 区 分 河道 内 部 冰 情 水 情 信息 ， 
但 该 指数 可 结合 历史 高 清 影像 作为 分 类 化 流程 实 
验 的 第 一 步 , 即 提取 河道 边界 ,排除 研究 对 象 外 的 
无 关 地 物 , 使 真实 河 冰 、 水 体 像 元 尽 最 大 限度 得 以 
保留 ”。 

黄河 内 蒙古 段 河 冰 信 息 提取 与 前 人 大 尺度 海 
冰 ( 湖 冰 ) 提 取 的 研究 相 比 , 既 具 统一 性 ,也 具 差 异 
性 。 例 如 , 海 勃 湾 库 区 是 典型 的 库 区 型 河 段 ,水 流 
流速 相对 较 小 , 河 冰 信息 提取 方法 同 海 冰 ( 湖 冰 ) 提 
取 大 体 相 同 ; 巴 并 高 勒 河 段 则 为 典型 的 游荡 型 河 
段 , 河 段 淤积 严重 , 河 冰 漫 滩 现 象 严重 , 河 冰 信 息 提 
取 受 不 确定 因素 增加 。5 种 遥感 指数 模型 在 黄河 内 
蒙古 段 冰 情 信 息 提取 精度 对 阔 值 的 响应 不 尽 相 同 ， 
E £L UR Bx sc 5 K RY ELTE Te BURG BEXT B (EE UII Jr e 
为 敏感 , 即 阀 值 的 少量 波动 就 会 对 提取 精度 产生 较 

影响 ;NDSI、MNDSI 以 及 NDWI 在 流程 化 实验 中 

提取 精度 对 阔 值 响应 适中 ,但 精度 略 小 于 近 红 外 波 
段 反射 率 阔 值 法 ;NDUWI 在 多 时 相 、 多 河 段 研究 中 
表现 出 较 高 提取 精度 ,同时 发 现 其 河 冰 水 体 最 优 
区 分 阔 值 大 体 分 布 于 姜 值 中 值 附 近 。 同 时 发 现 , 四 
岸 侧 模型 指数 值 略 大 于 上 四 岸 侧 ,这 可 能 与 弯 道 处 离 
心力 的 存在 导致 止 岸 侧 堆 冰 增加 有 关 , 故 河道 四 岸 
Ai ER EE E S p 779, 

黄河 冰 情 遥感 监测 不 仅 受 数 据 质量 FAR AR TIE 
以 及 光照 等 的 影响 ,并 且 也 受 水 位 ,流量 以 及 上 游 
水 库 防 凌 调 度 的 影响 ””。 例 如 :气温 越 高 , 河 冰 分 
布 范 围 越 小 , 稳 封 期 清 沟 所 占 比 例 越 大 ,已 发 生 凑 
汛 灾害 ;水 位 越 高 ,河道 与 大 气 表面 接触 面积 就 越 
大 , 河 冰 冻结 所 需要 的 热量 条 件 也 随 之 增 大 , 河 冰 
分 布 范 围 增 大 , 河 冰 漫 堤 风险 随 之 增加 等 。 
此 ,在 未 来 研究 中 ,如 何在 模型 中 结合 气象 条 件 .水 
利 条 件 等 数据 以 及 自动 化 提取 冰 水 像 元 ,同时 结合 
水 文 站 野外 观测 资料 做 对 比分 析 , 是 一 个 值得 考虑 
的 方向 ,这 可 能 会 进一步 提高 凌 汛 期 黄河 河 冰 水 
体 情 监测 的 时 效 性 和 准确 性 。 


4 结论 


黄河 内 蒙古 段 河 势 复 杂 、 地 形 地 貌 多 样 , 基 于 
遥感 影像 模型 监测 河 冰 、 水 体 情 分 布 不 易 获 取 且 各 
指数 模型 的 适用 性 尚未 有 人 做 出 评 佑 。 针 对 上 述 
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问题 ,本文 运用 L8 影 像 数据 ,通过 设计 实验 方案 提 
取 河 道 边界 ,利用 近 红 外 波段 反射 率 阔 值 NDWI、 
NDSI、MNDSI 及 NDUWI5 种 遥感 模型 对 研究 区 河 
UK 水体 分 布 进行 监测 ,进而 分 析 各 指数 的 阔 值 稳 
定性 及 不 确定 性 ,得 出 以 下 结论 

(1) 基于 NDSI 与 Coogle Earth 历史 高 清 影像 结 
合 的 二 次 排除 无 关 地 物 的 实验 方案 在 河 冰 \ 水 体 情 
监测 时 能 较 好 的 去 除 周围 无 关 地 物 。 

(2) NDUWI 在 黄河 内 蒙古 段 影 像 中 的 分 类 精 
度 均 达 到 90.00% 及 0.90 以 上 ,总 体 分 类 精度 较 高 ， 
适用 性 较 强 , 且 NDUWI 在 5 期 遥感 影像 中 的 河 冰 、 
水 体 最 优 区 分 阀 值 大 体 分 布 于 阔 值 中 值 附近 ;但 
NDUWI 无 法 将 土壤 与 河 冰 进行 有 效 区 分 ,使 用 前 应 
先 提取 河道 矢量 边界 ,排除 无 关 地 物 。 
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Ice-water classification in Inner Mongolia reach of the 
Yellow River based on remote sensing images 


ZHAI Yongguang, ZHANG Xin' JI Honglan, MOU Xianyou, ZHANG Baosen 
. College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, Hohhot , Inner mongolia, 
1. College of Water C y and Civil Engi ing, I Mongolia Agricultural University, Hohhot 010018, I goli 
China; 2. Yellow River Institute of Hydraulic Research, Zhengzhou 450003, Henan, China) 


Abstract: Timely delivery of detailed information on the spatial distribution of river ice during ice-flood season 
is highly valuable for predicting, and improving communication on, ice disaster. Remote sensing technology 
provides a key method for obtaining the spatial distribution of river ice. However, the large amount of suspended 
sediment in the Yellow River represents a challenge to high-precision discrimination between ice and water based 
on remote sensing technology. Taking the Inner Mongolia reach of the Yellow River as an example, this study 
compares and evaluates the performance of five indices in the classification of river ice and water: near-infrared 
reflectance; normalized difference water index (NDWI); normalized difference snow index (NDSI); improved 
normalized snow index (MNDSI); and normalized difference unfrozen water index (NDUWI). The overall 
classification accuracy and Kappa coefficient (a measure of reliability) were calculated for each index, and the 
threshold stability of each index was analyzed. The results show that NDUWI achieves the highest accuracy and 
reliability (Kappa coefficient) in each studied subregion. The overall classification accuracy and Kappa 
coefficient of NDUWI are more than 90.0% and 0.90, respectively, and the optimal discrimination threshold 
between river ice and ice-free water is close to the median value. These results can provide a basis for the selection 
of ice monitoring methods and optimization of ice-blasting locations on the Yellow River during ice-flood season. 
Key words: river ice; remote sensing index; threshold stability; uncertainty; Landsat 8; Inner Mongolia reach 
of the Yellow River 


